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Abstract. The hypergiant IRC +10420 is a unique object for 
the study of stellar evolution since it is the only object that is be- 
lieved to be witnessed in its rapid transition from the red super- 
giant stage to the Wolf-Rayet phase. Its effective temperature 
has increased by 1 000-2000 K within only 20 yr We present 
the first speckle observations of IRC +10 420 with 73 mas reso- 
lution. A diffraction-limited 2.11 /im image was reconstructed 
from 6 m telescope speckle data using the bispectrum speckle- 
interferometry method. The visibility function shows that the 
dust shell contributes ^ 40% to the total flux and the unre- 
solved central object ~ 60%. 

Radiative transfer calculations have been performed to 
model both the spectral energy distribution and visibility func- 
tion. The grain sizes, a, were found to be in accordance with 
a standard distribution function, n{a) ^ a^"^^, with a ranging 
between aj^m = 0.005 /im and Omax = 0.45 /^m. The observed 
dust shell properties cannot be fitted by single-shell models 
but seem to require multiple components. At a certain distance 
we considered an enhancement over the assumed l/r^ density 
distribution. The best model for both SED and visibility was 
found for a dust shefl with a dust temperature of lOOOK at its 
inner radius of 69 i?*. At a distance of 308 i?* the density was 
enhanced by a factor of 40 and and its density exponent was 
changed from a; = 2 to a; = 1.7. The shell's intensity distribu- 
tion was found to be ring-like. The ring diameter is equal to the 
inner diameter of the hot shell ('^69 mas). The diameter of the 
central star is ^ 1 mas. The assumption of a hotter inner shell 
of 1200 K gives fits of almost comparable quality but decreases 
the spatial extension of both shells' inner boundaries by ~ 30% 
(with X = 1.5 in the outer shell). The two-component model 
can be interpreted in terms of a termination of an enhanced 
mass-loss phase roughly 60 to 90 yr (for d — 5kpc) ago. The 
bolometric flux, i^boi, is 8.17 • 10^^° Wm^^ corresponding to 
a central-star luminosity of L/Lq — 25 462 • (d/kpc)^. 

Key words: Techniques: image processing — Circumstellar 
matter — Stars: individual: IRC +10420 — Stars: mass-loss 
— Stars: supergiants — Stars: Wolf-Rayet 



1. Introduction 

The star IRC +10420 (= V 1302 Aql = IRAS 19244+1115) is 
an outstanding object for the study of stellar evolution. Its spec- 
tral type changed from F8 1+ in 1973 (Humphreys et al. 1973) 
to mid-A today (Oudmaijer et al. 1996 Klochkova et al. 1997) 
corresponding to an effective temperature increase of 1000- 
2000 K within only 20 yr. It is one of the brightest IRAS ob- 
jects and one of the warmest stellar OH maser sources known 
(Giguere et al. |1976|, see also Mutel et al. |1979|, Diamond et 



al. 1983 Bowers 1984, Nedoluha & Bowers 1992). Ammonia 



emission has been reported by McLaren & Betz (1980) and 
Menten & Alcolea (1995). Large mass-loss rates, typically of 
the order of several 10"'' Mo/yr (Knapp & Morris |1985[, Oud- 



maijer et al. 1996 1 were determined by CO observations. Two 



evolutionary scenarios have been suggested for IRC +10420: 
It is either a post-AGB (AGB: Asymptotic Giant Branch) star 
evolvi ng thr ough the proto -plane tary nebula stage ( e.g. Fi x & 
Cobb |l987t Hrivnak et al. |l989t Bowers & Knapp |l989| ), or 
it is a massive hypergiant evolving from the RSG (Red Super- 
giant Branch) branch towards the Wolf-Rayet phase (e.g. Mutel 
et al. |l979t Nedol uha & Bowers |l992i Jo nes et al. |l993i Oud- 
maijer et al. 199f , Klochkova et al. 1997 ). However, due to its 
distance (rf = 3-5 kpc), large wind velocity (40 km/s) and photo- 
metric history, IRC +10420 is most likely a luminous massive 
star (see Jones et al. 1993 and Oudmaijer et al. 1996), there- 



fore being the only massive object observed until now in its 
transition to the Wolf-Rayet phase. The structure of the circum- 
stellar environment of IRC + 10420 appears to be very complex 
(Humphreys et al. 1997), and scenarios proposed to explain the 
observed spectral features of IRC + 10420 include a rotating 
equatorial disk (Jones et al. 1993), bipolar outflows (Oudmaijer 
et al. 1994), and the infall of circumstellar material (Oudmaijer 
1998). 

Previous infrared speckle and coronographic observations 



were reported by Dyck et al. (1984), Ridgway et al. (1986), 



Cobb & Fix (1987), Christou et al. (199C) and Kastner & 



Weintraub (1995 1. In this paper we present diffraction-limited 



73 mas bispectrum speckle-interferometry observations of the 
dust shell of IRC + 10 420 as well as radiative transfer calcula- 
tions to model its spectral energy distribution and visibility. 
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2. Observations and data reduction 

The IRC +10420 speckle interferograms were obtained with 
the Russian 6 m telescope at the Special Astrophysical Ob- 
servatory on June 13 and 14, 1998. The speckle data were 
recorded with our NICMOS-3 speckle camera (HgCdTe array, 
256^ pixels, frame rate 2 frames/s) through an interference fil- 
ter with a centre wavelength of 2.11 /j,m and a bandwidth of 
0.19 /im. Speckle interferograms of the unresolved star HIP 
95447 were taken for the compensation of the speckle inter- 
ferometry transfer function. The observational parameters were 
as follows: exposure time/frame 50 ms; number of frames 8400 
(5200 of IRC + 1 420 and 3200 of HIP 95447); 2. 1 1 ^m seeing 
(FWHM) -I'.'O; field of view 7'.'8x7'.'8; pixel size 30.5 mas. A 
diffraction-limited image of IRC +10420 with 73 mas resolu- 
tion was reconstructed from the speckle interferograms using 
the bispectrum speckle-interferometry method (Weigelt 1977] , 
Lohmann et al. 1983, Hofmann & Weigelt 1986). The bispec- 
trum of each frame consisted of ^37 million elements. The 
modulus of the object Fourier transform (visibility) was de- 
termined with the speckle interferometry method (Labeyrie 



197C). 



It is noteworthy that 2.11 /im filters also serve to im- 
age hydrogen emission as, for instance, H2 (2.125/im) or 
Br 7 (2.166/im) emission. Accordingly, it is possible that one 
might look at hydrogen emission rather than at the dust emis- 
sion of a circumstellar shell. The low-resolution spectr um o f 
IRC +10 420 published in the atlas of Hanson et al. ( |1996{ ) 
shows a Br 7 line in emission as the most prominent feature for 
the wavelength range considered here. Oudmaijer et al. (1994) 
carried out high-resolution infrared spectroscopy and found an 
equivalent width of 1 .2 A for the Br 7 emission line. This is 
only 0.06% of the bandwidth of our interference filter and con- 
sequently negligible. 

Figure [l| shows the reconstructed 2.11 /im visibility func- 
tion of IRC +10420. There is only marginal evidence for an 
elliptical visibility shape (position angle of the long axis ^ 
130° ± 20°, axis ratio - 1.0 to 1.1). The visibility 0.6 at fre- 
quencies > 4 cycles/arcsec shows that the stellar contribution 
to the total flux is -^60% and the dust shell contribution is 
^40%. In order to compare our results with speckle obser- 
vations of other groups we determined the GauB fit FWHM 
diameter of the dust shell to be dpwHM — (219 ± 30) mas. 
By comparison, Christou et al. (1990) found for 3.8 m tele- 
scope K-band data a dust-shell flux contribution of '^50% 
and rfpwHM — 216 mas. However, as will be shown later, 
a ring-like intensity distribution appears to be much better 
suited than the assumption of a Gaussian distribution whose 
corresponding FWHM diameter fit may give misleading sizes 
(see Sect. 3.4.5). Fig. displays the azimuthally averaged 



diffraction-limited images of IRC +10420 and the unresolved 
star HIP 95447. 



3. Dust shell models 
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Fig. 1. Left: Two-dimensional 2.11 /^m visibility function of 
IRC + 10420 shown up to the diffraction limit (see right panel). 
The dark central structure shows that the central object is sur- 
rounded by a dust shell. Right: Azimuthally averaged 2.11 /im 
visibility of IRC +10420 with error bars for selected frequen- 
cies. This visibility function consists of a constant plateau (vis- 
ibility ^ 0.6) caused by the unresolved central object and a 
triangle-shaped low-frequency function caused by the faint ex- 
tended nebula. 



200 



Fig. 2. Azimuthally 
averaged radial plots 
of the reconstructed 
diffraction-limited 
2.11/im-images of 
IRC + 10420 (solid 
line) and HIP 95447 
(dashed line). 




mage Space (mas) 



3.1. Spectral energy distribution 

The spectral energy distrubion (SED) of IRC +10 420 with 
9.7 and 18 /im silicate emission features is shown in Fig. ^ 
It c orresp onds to the '1992' data set used by Oudmaijer et 
al. ( 1996 ) and combines VRI (October 1991), near-infrared 
(March and April 1992) and Kuiper Airborne Observatory pho- 
tometry (June 1991) of Jones et al. ( |l993| ) with the IRAS mea- 
surements and 1.3 mm data from Walmsley et al. (1991). Ad- 
ditionally, we included the data of Craine et al. (1976) for 
A < 0.55 /im. In contrast to the near-infrared, the optical mag- 
nitudes have remained constant during the last twenty years 
within a tolerance of w 0'."1. 

IRC + 10420 is highly reddened due to an extinction of 
7™ by the interstellar medium and the circumstellar 
shell. From polarization studies Craine et al. (1976) estimated 
an interstellar extinction of Ay « 6™. Jones et al. (1993) de- 
rived from their polarization data Ay « 6'" to 7™. Based on 
the strength of the diffuse interstellar bands Oudmaijer ( 1998| ) 
inferred E{B - V) = r."4 ± ff^S for the interstellar contri- 
bution compared to a total of E{B — V) — 2™4. We w ill use 
an interstellar Ay of 5™ as in Oudmaijer et al. ( 1996| ). This 



/I total 
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interstellar reddening was taken into account by adopting the 
method of Savage & Mathis (|1979|) with Ay = 3.1E{B - V). 



3.2. The radiative transfer code 

In order to model both the observed SED and 2.11 ^va visibil- 
ity, we performed radiative transfer calculations for dust shells 
assuming spherical sym metry. We used the code DUSTY de- 
veloped by Ivezic et al. ( 1997 ), which solves the spherical ra- 
diative transfer problem utilizing the self-similarity and scaling 
behaviou r of IR emission from radiatively heated dust (Ivezic 
& Elitzur 1997). To solve the radiative transfer problem includ- 
ing absorption, emission and scattering several properties of the 
central source and its surrounding envelope are required, viz. 
(i) the spectral shape of the central source's radiation; (ii) the 
dust properties, i.e. the envelope's chemical composition and 
grain size distribution as well as the dust temperature at the in- 
ner boundary; (iii) the relative thickness of the envelope, i.e. the 
ratio of outer to inner shell radius, and the density distribution; 
and (iv) the total optical depth at a given reference wavelength. 
The code has been expanded for the calculation of synthetic 
visibilities as described by Gauger et al. (1999). 



3.3. Single-shell models 

We calculated various models considering the following pa- 
rameters within the radiative transfer calculations: SED and 
visibility were modelled for Toff = 7000 to 9000 K, black bod- 
ies and Kurucz (1992) model atmospheres a s cen tral sources 
of radiation, different silicates (Draine & Lee 1984 , Ossenkopf 
et al. 1992, David & Pegourie 1995), single-sized grains with 



a = 0.01 to 0.6 /im and grain size distributions according to 
Mathis et al. ( |1977[ hereafter MRN), i.e. n{a) ^ a^^-^, with 
0.005 /^m < a < (0.20 to 0.60) ^m. We used a l/r^ density 
distribution and a shell thickness Fout = ?'out/?'i of 10'^ to 10^ 
with Tout and Ti being the outer and inner radius of the shell, 
respectively. Then, the remaining fit parameters are the dust 
temperature, Ti, which determines the radius of the shell's in- 
ner boundary, ri, and the optical depth, r, at a given reference 
wavelength, Arof- We refer to Arof = 0.55 /im. Models were 
calculated for dust temperatures between 400 and lOOOK and 
optical depths between 1 and 12. Significantly larger values for 
T lead to silicate features in absorption. 

Fig. I shows the SED calculated for Teff=7000K, 
yoyt=10^, Draine & Lee (1984) siHcates, MRN grain size dis- 
tribution (flmax — 0.2 /im) and different dust temperatures. It 
illustrates that the long-wavelength range is sufficiently well 
fitted for cool dust with Ti = 400 K, optical wavelengths and 
silicate features require t ^ 5. The inner radius of the dust shell 
is at ri = 447 i?* (i?* : stellar radius), the equilibrium temper- 
ature at the outer boundary amounts to Tout — 22 K. However, 
the fit fails in the near-infrared underestimating the flux be- 
tween 2 and 5 /im. Instead this part of the SED seems to require 
much hotter dust of Ti 1> 800K(ri < 145 R ^ Tpu t = 32 K). 
This confirms the findings of Oudmaijer et al. ( 1996| ) who con- 
ducted radiative transfer calculations in the small particle limit. 




X [urn] 

Fig. 3. Model SED for Tcff = 7000 K, ro.ssMm = 5 and dif- 
ferent dust temperatures Ti . The lower panel shows the silicate 
features. The calculations are based on a black body, Draine & 
Lee (1984) silicates, and an MRN grain size distribution with 
flmax = 0.2 /im. The symbols {+) refer to the observations (see 
text) corrected for interstellar extinction of A^ — 5"\ 




Fig. 4. Model visibihty function for Tcff = 7000 K, To.ss^m = 
5, Ti = 600 K and different maximum grain sizes in the MRN 
grain size distribution (omax — 0.2, 0.4, 0.45, 0.5, 0.55 and 
0.6 /im from top to bottom). The calculations are based on a 
black body and Draine & Lee (1984) silicates. 

where scattering is negligible. They introduced a cool (400 K) 
and a hot (lOOOK) shell to achieve an overall fit. We found 
this behaviour of single-shell SEDs to be almost independent 
of various input parameters. Increasing l^ut from 10'^ to 10^ 
leads to somewhat higher fluxes, but only for A > 100/im. 
The equilibrium temperature at the outer boundary decreases 
by a factor of two if the shell's thickness is increased by one 
order of magnitude. Larger Tcff gives slightly less flux in the 
near-infrared, larger wavelengths (A > 10/im) are almost unaf- 
fected. The Draine & Lee (1984) and David & Pegourie (1994) 
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silicates give almost identical results, the optical constants of 
Ossenkopf et al. (1992) lead to a larger 9.7 /im/18 /im flux ra- 
tio for the silicate features, to somewhat higher fluxes between 
2 and 10 /im and to a somewhat flatter slope of the SED at 
short wavelengths. However, the need for two dust compo- 
nents still exists. Calculations with different grain sizes show 
that single-sized grains larger than 0.2 /im are not suitable for 
IRC H- 10 420. The silicate features are worse fitted and, in par- 
ticular, a significant flux deficit appears in the optical and near- 
infrared. The variation of the maximum grain size in the MRN 
distribution leads to much smaller differences due to the steep 
decrease of the grain number density with grain size. 

The 2.11 ^m visibility is very sensitive against scattering, 
thus depending strongly on the assumed grain sizes (see Groe- 
newegen 1997 ) as demonstrated in Fig. ^. For a given set of 
parameters both inclination and curvature of the visibility are 
mainly given by the optical depth, r, and the grain size, a. Since 
r is fixed to small values due to the emission profiles, a can 
be determined. The dust temperature must be varied simulta- 
neously since an increase of Ti leads to a steeper declining 
visibility. Our calculations show that the visibility is best fit- 
ted for an intermediate Ti = 600 K in contrast to the SED. 
Either single-sized grains with a ^ 0.4 /im (which, however, 
are ruled out by the SED) or MRN grain size distributions 
with flinax '^ 0.45 to 0.5 /im are appropriate. This result still 
depends on the kind of silicates considered, i.e. on the opti- 
cal constants. For instance, if we take the 'warm silicates' of 
Ossenkopf et al. (1992), we get somewhat smaller particles 
(by ~ 0.1 /im, i.e. a ^ 0.3 /im for single-sized grains and 
Omax ~ 0.35 /im for a grain distribution, resp.). The differ- 
ences to the corresponding 'cold siUcates' or to the data from 
David & Pegourie (1994) are found to be smaller. The fits to the 
SED are of comparable quality. We chose Draine & Lee ( 1984 ) 
silicates with Omin = 0.005 /im and Omax — 0.45 /im. 

3.4. Multiple dust-shell compenents 



3.4.1. Two component shells 

Since we failed to model the SED with the assumptions made 
so far, we introduced a two-component shell as Oudmaijer et 
al. (1996). For that purpose, we assume that IRCh-10420 had 
passed through a superwind phase in its history as can be ex- 
pected from its evolutionary status (see Schaller et al. 1992, 
Garcia-Segura et al. 1996). This is in line with the conclusions 
drawn from the Oudmaijer et al. (1996) model and recent in- 
terpretations of HST data (Humphreys et al. 1997). A previous 
superwind phase leads to changes in the density distribution, 
i.e. there is a region in the dust shell which shows a density en- 
hancement over the normal 1/r^ distribution. The radial den- 
sity distribution may also change within this superwind shell. 
For more details, see Suh & Jones ( 1997[ ). Since dust formation 
operates on very short timescales in OH/IR stars, we assume a 
constant outflow velocity for most of the superwind phase and 
thus a 1/r^ density distribution. For simplicity, we consider 
only single jumps with enhancement factors, or amplitudes, A 
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Fig. 5. Relative density dis- 
tribution for a superwind at 
Y = r/ri = 4.5 with an 
amplitude of A = 10. 




Fig. 6. SED (top) and visibility (bottom) for a superwind 
model with Y = r/ri = 4.5 and different amplitudes. Model 
parameters are: black body. Toff = 7000 K, Ti = 1000 K, 
^0.55^111 — 7.0, Draine & Lee (1984) silicates, Mathis et 
al. (1977) grain size distribution with flmax ~ 0.45 /im, and 
i^out = 10"*. The symbols refer to the observations (see text) 
corrected for interstellar extinction of Ay = 5"\ 



at radii Y = r/ri in the relative density distribution as demon 
strated in Fig. pi 



Concerning the grains we stay with Draine & Lee ( |1984| ) 
silicates and an MRN grain size distribution with Omin = 
0.005 /im and amax = 0.45 /im as in the case of the single 
shell models. The influence of different grain-size distributions 
will be discussed later 

We calculated a grid of models for Ti = 1000 K with su- 
perwinds at F = 2.5 to 8.5 with amplitudes A ranging from 10 
to 80. Due to the introduced density discontinuity the flux con- 
servation has to be controlled carefully, in particular at larger 
optical depths and amplitudes. SED and visibility behave con- 
trarily concerning the adjustment of the superwind: The SED 
requires sufficiently large distances, Y <; 4.5, and moderate 
amplitudes, A < 20 to 40, in particular for the flux between 
2 and 10 /im and for A > 20/im. A good fit was found for 
y = 6.5 and A = 20 corresponding to ri = 81 i?*. Note that 
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the bolometric flux at the inner dust-shell radius (and there- 
fore ri/r,) is fully determined by the solution of the radia- 
tive transfer problem even though the overall luminosity is not 
(Ivezic & Elitzur 1997). The dust temperature at the density 
enhancement (r2 = 527 i?*) has dropped to 322 K. This agrees 
well with the model of Oudmaijer et al. (1996). The visibil- 
ity, however, behaves differently. In order to reproduce the un- 
resolved component (the plateau) large amplitudes, yl > 40 
to 80, are required. On the other hand, the slope at low spa- 
tial frequencies is best reproduced for a close superwind shell, 
Y < 4.5 (at this distance independent on A). The best model 
found for both SED and visibility is that with Y = 4.5 and 
A = 40 as shown in Fig. ^. It corresponds to ri = 71 i?* 
and r2 = 320 i?* (with T2 ^ 475 K), i.e. to angular diam- 
eters of 01 = 71 mas and 82 = 321 mas. The angular di- 
ameters depend on the model's bolometric flux, i^boh which 
is 8.17 • 10"^" Wm^^. Accordingly, the central star has a lu- 
minosity of L/Lq ~ 25 462 • (d/kpc)^ and an angular di- 
ameter of e* = 1.74 ■ 10^ ^/F\~JT^ - Imas. Assum- 
ing a constant outflow velocity of w = 40km/s, the expan- 
sion ages of the two components are ti/yr = 4.2 • (d/kpc) 
and t2/yr = 18.9 • (d/kpc). With a dust-to-gas ratio of 0.005 
and a specific dust density of 3 gcm~^ the mass-loss rates of 
the components are Mi = 1.4 • lO~^M0/yr • (d/kpc) and 
A/2 = 5.5 • IO-^Mq/yv ■ (d/kpc). 

Fig. n\ shows the fractional contributions of the direct stel- 
lar radiation, the scattered radiation and the dust emission to 
the total emerging flux. The stellar contribution has its max- 
imum at 2.2 /im where it contributes 60.4% to the total flux 
in accordance with the observed visibility plateau of 0.6. At 
this wavelength scattered radiation and dust emission amount 
to 25.6% and 14% of the total flux, respectively. Accordingly, 
64.6% of the 2.11 /im dust-shell emission is due to scattered 
stellar light and 35.3% due to direct thermal emission from 
dust. For A Si 1 fim the flux is determined by scattered radiation 
whereas for A ^ 10 /im dust emission dominates completely. 



3.4.2. Influence of the grain-size distribution 

As in the case of the single-shell models we also studied other 
grain size distributions. The MRN distribution derived for the 
interstellar medium gives a continuous decrease of the num- 
ber density with increasing grain sizes. On the other hand, the 
distribution of grains in dust-shells of evolved stars rather ap- 
pears to be peaked at a dominant size (e.g. Kriiger & Sedlmayr 
1997, Winters et al. 1997). It is noteworthy that even in the case 
of a sharply peaked size distribution the few larger particles 
can contribute significantly to the absorption and scattering co- 
efficients (see Winters et al. 1997). Accordingly, the 2.11 /im 
visibility reacts sensitively if some larger particles are added 
whereas the SED does not, as demonstrated in the previous 
section. In order to study the influence of different grain size 
distributions on the two-component model we calculated grids 
of models with n{a) ^ a'* for different exponents (q = —3.0 to 
—5.5) and lower and upper cut-offs (amin = 0.005 to 0.05 /im 
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Fig. 7. Fractional contributions of the emerging stellar radiation 
as well as of the scattered radiation and of the dust emission to 
the total flux as a function of the wavelength for a superwind 
model with Y — r/ri = 4.5 and A — 40. Model parameters 
are: black body. Toff = 7000 K, Ti = 1000 K, ro.ss^m = 7.0, 
Draine & Lee (1984) silicates, Mathis et al. (1977) grain size 



distribution with a„ 



0.45 /im. 



and flinax ~ 0.1 to 0.8 /im). Additionally we considered single- 
sized grains (a = 0.1 to 0.8 /im). 

Concerning the visibility, a larger (smaller) negative ex- 
ponent in the distribution function can, in principle, be com- 
pensated by increasing (decreasing) the maximum grain size. 
For instance, q = —4.0 requires Omax = 0.55 /im to fit the 
2.11 /im visibility. On the other hand, if the distribution be- 
comes too narrow, the SED cannot be fitted any longer since 
the 9.7 /im silicate feature turns into absorption. A distribution 
with q = —3.8 and Omax — 0.50 /im best reproduces the flux- 
peak ratio of the silicate features. 

For a given exponent in the grain-size distribution function 
of q = —3.5 we arrive at the same maximum grain size as in 
the case of the one-component model, viz. 0.45 /im, in order to 
yield a fit for both the SED and the visibility (see Fig. M). This 
is due to the fact that larger particles increase the curvature of 
the visibility curve at low spatial frequencies whereas the high- 
frequency tail (the plateau) is found at lower visibility values. 
On the other hand, the inclusion of some larger particles does 
not change the shape of the SED as discussed above. 

If sufficiently small, the lower cut-off grain size can be 
changed moderately (within a factor of two) without any sig- 
nificant change for SED and visibility. If amin exceeds, say, 
0.05/im, the fits of the observations begin to become worse. 
For instance, the curvature of the visibility at low spatial fre- 
quencies and the flux-peak ratio of the silicate features are then 
overestimated. 

Finally, we repeated the calculations under the assumption 
of single-sized grains. In order to model the visibility a grain 
size a close to 0.3 /im is required as shown in Fig. pi In con- 
trast, the reproduction of the relative strengths of the silicate 
features seems to require smaller grains, viz. close to 0.1 /im. 
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Fig. 8. SED (top) and visibility (bottom) for a superwind 
model with Y = r/ri = 4.5 and A = 40 calculated for Mathis 
et al. (1977) grain size distributions with Omax = 0.20, 0.4, 0.5 
and 0.6 /im. Model parameters are: black body, Teff = 7000 K, 
Ti = 1000 K, ro.55^m = 7.3, Draine & Lee (1984) siHcates, 
and Yout = 10"*. The symbols refer to the observations (see 
text) corrected for interstellar extinction of A^ — 5™. 



Fig. 9. SED (top) and visibility (bottom) for a superwind 
model with Y = r/ri = 4.5 and A = 40 for single-sized 
grains with a = 0.1, 0.2 and 0.3 /im. Model parameters are: 
black body, Teff = 7000 K, Ti = 1000 K, ro.ssMm = 7.3, 
Draine & Lee (1984) silicates, and l^ut = 10^. The symbols 
refer to the observations (see text) corrected for interstellar ex- 
tinction of ^v = 5™. 



Consequently, for the modelling of IRC H- 10 420 a grain size 
distribution appears to be much better suited than single-sized 
grains. 

3.4.3. Influence of the density distribution 

Inspection of the best fits derived so far reveals that there are 
stiU some shortcomings of the models. First, although being 
within the observational error bars, the model visibilities al- 
ways show a larger curvature at low spatial frequencies. This 
seems to be almost independent of the chosen grain-size dis- 
tribution. Second, the flux beyond 20 /im is somewhat too low. 
This may be due to our choice of a l/r^ density distribution for 
both shells. We recalculated the model grid for different 1/r^ 
density distributions for both shells with x ranging between 1 
and 4. A flatter distribution in the outer shell increases the flux 
in the long-wavelength range as required but leads also to a 
drop of the flux in the near-infrared. The plateau in the visi- 
bility curve remains unaffected but the curvature at low spatial 
frequencies is increased. To take advantage of the better far- 
infrared properties of cool shells with flatter density distribu- 
tions, but to counteract their disadvantage in the near-infrared 
and at low spatial frequencies, the density distribution of the in- 
ner shell also has to be changed. It should be somewhat steeper 
than the normal 1/r^ distribution. Then the near-infrared flux 
is raised and the visibility shows a smaller curvature in the 



low-frequency range. It should be noted that the curvature is 
most affected for superwinds of low amplitudes. However, the 
steeper density decrease in the inner shell leads to increasingly 
low visibility values in the high frequency range. Since this has 
to be compensated by an increase of the superwind amplitude 
the advantages of the steeper distribution are almost cancelled. 

Thus, we can stay with a 1/r'^ density distrubution in the 
inner shell and moderate superwind amplitudes (A ^ 40). The 
then best suited models we found are those with superwinds at 
Y = 4.5 and a 1/r^-^ distribution in the outer shell. The corre- 
sponding SED and visibility are shown in Fig. HQ for different 
superwind amplitudes. We note again that the quality of the fits 
is in particular determined by the outer shell, whereas the in- 
ner shell's exponent is less constrained. A 1/r^ distribution in 
the inner shell and large superwind amplitudes (A <; 80) give 
similar results. 

The radii of the inner and outer shell are ri — 69 R^ and 
r2 = 308 i?* (with T2 ~ 483 K), resp., corresponding to angu- 
lar diameters of 81 = 69 mas and 82 = 311 mas. Adopting the 
same assumptions for outflow velocity, dust-to-gas ratio and 
specific dust density as in the previous section, the expansion 
ages are ti/yr — 4.1 ■ (d/kpc) and t2/yr — 18.4 • (d/kpc), 
for the mass-loss rate of the inner component one gets Mi = 
1.35 ■ 10^^ Mq/jt ■ (d/kpc). In the outer component either 
the outflow velocity has increased or the mass-loss rate has 
decreased with time due to the more shallow density distri- 
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Fig. 10. SED (top), silicate features (middle) and visibility 
(bottom) for a superwind model with Y ~ r/ri — 4.5 and 
different amplitudes. The inner shell obeys a 1/r^ density dis- 
tribution, the outer shell a l/r^^ density distribution. Model 
parameters are: black body. Toff = 7000 K, Ti = 1000 K, 
To.55fim = 7.0, Draine & Lee (1984) silicates, Mathis et 
al. (1977) grain size distribution with a,nax — 0.45 /im, and 
^out = 10"*. The symbols refer to the observations (see text) 
corrected for interstellar extinction of A^ — 5™. 



bution. Provided the outflow velocity has kept constant, the 
mass-loss rate at the end of the superwind phase, 92 yr ago, was 
M2 — 5.4-10^* AIq/yy- (d/kpc), and, for instance, amounted 
to M2 = 8.0 • 10-4 Mq/jt ■ (d/kpc) 200 yr ago. 

Since the flatter density distribution provides a better fit for 
the long-wavelength range of the SED, while the visibility is 
equally well fitted compared to the standard density distribu- 
tion, it is superior to the model of Sect. 3.4.1 Fig. O gives 



the fractional flux contributions (stellar, dust, scattering) for the 
same model as shown in Fig. p^but with an l/r^-^ distribution 
in the outer shell. The various flux contributions at 2. 1 1 /im are 
very similar to those of the 1/r^ model: 62.2% stellar light, 
26.1% scattered radiation and 10.7% dust emission. Thus, the 
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Fig. 11. Fractional contributions of the emerging stellar radia- 
tion as well as of the scattered radiation and of the dust emis- 
sion to the total flux as a function of the wavelength for a super- 
wind model with Y — r/ri = 4.5, A — 40 and a l/r^-^ den- 
sity distribution in the outer shell. Model parameters are: black 
body, Teff = 7000 K, Ti = 1000 K, rcss^m = 7.0, Draine & 
Lee (1984) silicates, Mathis et al. (1977) grain size distribution 
with Cmax — 0.45 fim. 



total emission of the circumstellar shell is composed of 70.9% 
scattered stellar light and 29.1% direct thermal emission from 
dust. 



3.4.4. Influence of the dust temperature 

Finally, we studied the influence of the dust temperature at the 
inner boundary of the hot shell. For that purpose we recalcu- 
lated the previous model grids for dust temperatures of 800 
and 1200K. As akeady shown for the single-shell models, an 
increase of the temperature at the inner boundary increases the 
flux in the near-infrared and substantially lowers the flux in the 
long-wavelength range. On the other hand, the higher the tem- 
perature the less is the curvature of the visibility at low spatial 
frequencies, the plateau is only significantly affected for low- 
amplitude superwinds. The shape of SED and 2.11/im visi- 
bility for different dust temperatures at the hot shell's inner 
boundary for a given superwind is demonstrated in Fig. |l2[ 
A temperature less than 1000 K can be excluded in particu- 
lar due to the worse fit of the visibility for low frequencies. 
Instead, the 1200 K model gives a much better fit to the vis- 
ibility than previous ones. Fig. u2 refers to an amplitude of 
A = 40 in order to be directly comparable with the models 
shown before. We note that we get an even better fit assuming 
A = 80, which leaves the low-frequency -range unchanged but 
improves the agreement with the measured plateau. However, 
the improvement of the 2.11 /im visibility model due to a hotter 
inner shell with Ti = 1200 K is at the expense of a consider- 
able amplification of the flux deficit for A <; 20 /im in the SED. 
In order to compensate this effect we have had to assume a 
flatter density profile for the outer shell than in the case of the 
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Fig. 12. SED (top) and visibility (bottom) for a superwind 
model {Y = r/ri = 4.5 and A = 40) with different tem- 
peratures for the inner boundary of the hot shell. Model param- 
eters are: black body, T^s — 7000 K, To.55pm — 7.0, Draine & 
Lee (1984) silicates, Mathis et al. (1977) grain size distribution 
with flmax = 0.45 /im and Yout — 10"*. The symbols refer to 
the observations (see text) corrected for interstellar extinction 
of Av = 5'". 

Ti = 1000 K., viz. - l/ri-5 instead of - l/r^-'^. The corre- 
sponding curves are shown in Fig. [IJ. Again, increasing the 
far-infrared fluxes, as required to model the SED, leads to an 
increase of the 2.11 ^m visibility's curvature at low spatial fre- 
quencies giving somewhat worse fits for the visibility. We note 
that the peak-ratio of the silicate features is better matched with 
a lower dust temperature of Ti = 1000 K. 

The radii of the inner and outer shell are now considerably 
smaller than those of the previous models due to the higher 
temperature of the hot shell. The radiative transfer calculations 
give here n = 47 i?* and r2 = 210 i?* (with T2 - 594 K), 
resp., resulting in angular diameters of 0i =47 mas and 02 ~ 
212 mas. Accordingly, the expansion ages are ii/yr = 2.8 • 
(d/kpc) and t2/yr = 12.6 • (d/kpc), for the mass-loss rate 
of the inner component one gets Mi = 9.2 • 10^^ Mq/yi ■ 
(d/kpc). Provided the outflow velocity has kept constant, the 
mass-loss rate at end of the superwind phase, 63 yr ago, was 
M2 — 3.7-10^'* Afo/yr- (d/kpc), and, for instance, amounted 
to M2 = 5.4 • 10-4 Afo/yr • (d/kpc) 200 yr ago. 



3.4.5. Intensity distributions 

Fig. [redisplays the spatial distribution of the obtained normal- 
ized model intensity for the model shown in Fig. |l^ (Ti = 
1000 K, Y = 4.5, A = 40, l/r^ and l/r^'^ density distiibu- 



Fig. 13. SED (top) and visibility (bottom) for a superwind 
model with Y = r/ri = 4.5 and different amplitudes. The 
inner shell obeys a 1/r^ density distribution, the outer shell 
a 1/r^-^ density distribution. The temperature at the inner 
boundary of the hot shell is 1200 K. Model parameters are: 
black body, T^s = 7000 K, To.sg^m = 7.0, Draine & Lee 
(1984) silicates, Mathis et al. (1977) grain size distribution with 
flmax — 0.45 /im, and Yout — 10'*. The symbols refer to the 
observations (see text) corrected for interstellar extinction of 
Av = 5"\ 



tion in the inner and outer shell, resp.) The (unresolved) cen- 
tral peak belongs to the central star, and the two local intensity 
maxima to the loci of the inner rims of the two shells at 35 mas 
and 157 mas, resp. The 2.11 /im intensity shows a ring-like dis- 
tribution with a steep decline with increasing distance from the 
inner boundary of the circumstellar shell. Similarly shaped in- 
tensity distributions have also been found by Ivezic & Elitzur 
(1996) for optically thin shells. 

We recall that this intensity distribution is based on radia- 
tive transfer models taking into account both the SED and the 
2.11 /im visibility. Figure O shows the model visibilities for 
much higher spatial frequencies than covered by the present 
observations. The required baselines would correpond to ^ 22 
and 440m instead to 6m (13.6cycles/arcsec). Since the dust- 
shell's diameter is ^ 70 mas a plateau is only reached for spa- 
tial frequencies larger than, say, 15 cycles/arcsec depending 
on the strength of the superwind. The central star is resolved 
at spatial frequencies of ^ 1000 cycles/arcsec. At frequencies 
^ 2 cycles/arcsec the shape of the observed and the modelled 
visibility function is triangle-shaped, which is a consequence 
of the ring-like intensity distribution of the dust shell. 

Visibility observations are often characterized by fits 
with Gaussian intensity distributions. The resulting Gaussian 
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Fig. 14. Normalized intensity vs. angular displacement d for 
a superwind model with Y — r/ri = 4.5, A — 40 and 
a l/r^-^ density distribution in the outer shell. The (unre- 
solved) central peak belongs to the central star. The inner hot 
rim of the circumstellar shell has a radius of 35 mas, and the 
cool component is located at 155 mas. Both loci correspond 
to local intensity maxima. Model parameters are: black body, 
Teff = 7000 K, Ti = 1000 K, ro.sspm = 7.0, Draine & Lee 
(1984) silicates, Mathis et al. (1977) grain size distribution with 
Omax = 0.45 /im. 




A=10 Y=4.5 
A=20 Y=4.5 
A=40 Y=4.5 
A=80 Y=4.5 




400 600 

q [arcsec'^] 



1000 



Fig. 15. Model visibility up to 50 (top) and 1000 cycles/arcsec 
(bottom) for a superwind model with Y = r/ri = 4.5 and 
different amplitudes. The inner shell obeys a 1/r^ density dis- 
tribution, the outer shell a l/r^^ density distribution. Model 
parameters are: black body. Toff = 7000 K, Ti = 1000 K, 
To.55/jm = 7.0, Draine & Lee (1984) silicates, and Mathis et 
al. (1977) grain size distribution with Omax = 0.45 /xm. 



FWHM diameter is then assumed to give roughly the typical 
size of the dust shell. A GauB fit to the observed visibility 
would yield a FWHM dust-shell diameter of (219 ± 30) mas in 
agreement with the one given by Christou et al. (1990). How- 
ever, radiative transfer models show that a ring-like intensity 
distributions appears to be more appropriate than a Gaussian 
one for the dust shell of IRC H- 10 420. The distribution shows a 
limb-brightenend dust condensation zone and a ring diameter 
of 70 mas. 

4. Summary 

Radiative transfer calculations show that the near-infrared vis- 
ibility strongly constrains dust shell models since it is, e.g., 
a sensitive indicator of the grain size. Accordingly, high- 
resolution interferometry results provide essential ingredients 
for models of circumstellar dust-shells. Assuming spherical 
symmetry we carried out radiative transfer calculations for the 
hypergiant IRC H- 10 420 to model both its SED and 2.11 /im 
visibility. Since we failed to find good SED fits for single- 
component models, we improved our density distribution intro- 
ducing a second component with enhanced values at a certain 
distance. For different scaled distances Y = r/ri and density 
enhancements A of this cool component we considered differ- 
ent grain-size distributions n{a) ^ a', density distributions 
p ~ l/r^ within the shells, and temperatures Ti at the inner 
boundary of the hot shell. 



An MRN grain size distribution 71(a) 



,-3.5 



with 

0.005 /im < a < 0.45/im was found to be well suited for 
IRC H-10420. Larger negative exponents, i.e. a narrower distri- 
bution, can be accounted for by increasing the maximum grain 
size. For instance, n{a) ^ a^'^-^ requires Omax ^ 0.55 /im. 
However, the range of appropriate exponents seemed to be 
quite small and steeper declining distributions led to signifi- 
cantly worse fits. 

Assuming a l/r^ density distribution for both shells and 
Ti = 1000 K gives the best fit for y = 4.5 and A := 40 (Fig. |). 
This model can be improved by introducing a somewhat flatter 
density distribution, viz. ~ l/r^^, for the outer shell leading 
to a better match with the observed SED for A <; 20 /im. The 
quality of the visibility fit remains almost unchanged (Fig. nOb. 
Both models show a somewhat larger curvature of the visibility 
at low spatial frequencies. However, the deviations are within 
the observational uncertainties. The various flux contributions 
at 2.11/im are 62.2% stellar light, 26.1% scattered radiation 
and 10.7% dust emission. 

Alternatively one may increase the temperature at the inner 
boundary of the hot shell to Ti = 1200 K which gives some- 
what better matches to the near-infrared flux and lowers the 
low-frequency visibility curvature. To counteract the concomi- 
tant loss of flux in the far-infrared one has to assume a l/r^^ 
density distribution (Fig. |l3|). The fit to the silicate features is, 
however, somewhat worse than in the case of the Ti = 1000 K 
model. 

The intensity distribution was found to be ring-like. This 
appears to be typical for optically thin shells (here To.ss^m — 7, 
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''■2.11pm = 0.55; see also Ivezic & Elitzur 1996) showing limb- 
brightened dust-condensation zones. Accordingly, the interpre- 
tation of the observational data by FWHM GauB diameters may 
give misleading results. 

The two components can be interpreted as if IRC H- 10 420 
has suffered from much higher mass-loss rates in its recent 
past than today. For instance, the Ti = 1000 K model gives 
Ml = 1.4-10-5 Mo/yr-(d/kpc) and Af2 = 8.5-lO-'*Af0/yr- 
(d/kpc). The kinematic age of the outer component gives 
a corresponding timescale of ^ 100 yr (for d = 5kpc). If 
T2 = 1200 K both shells are located closer to the central star 
by approximately 30% leading to a correspondingly smaller 
timescale. The failure of constant mass-loss wind models to fit 
the SED agrees with the findings of Oudmaijeret al. (1996) and 
Humphreys et al. (1997). A previous high mass-loss episode is 
in line with the suspected post-RSG stage of IRC H- 10 420. 

Although, the present observations give only marginal evi- 
dence for deviation from spherical symmetry (if elliptical, po- 
sition angle of the long axis ^ 130° ± 20°, axis ratio ^ 1.0 to 
1.1), the hot shell may also be interpreted as a disk with a typ- 
ical diameter of approximately 50 mas. The presence of a ro- 
tating equatorial disk has been proposed by Jones et al. (1993), 
and Oudmaijer et al. (1996) interpreted their hot dust-shell as a 
disk as well. Provided the disk is not viewed pole-on, the cor- 
responding two-dimensional power spectra should be clearly 
elongated. It should be noted, however, that disks with an ex- 
tension of typically 50 mas can only be detected in the power 
spectra if they provide at least, say, 10% of the total flux. Oud- 
maijer (1995) discussed several models for the circumstellar 
shell of IRC H- 10 420 and found neither a bipolar nor a disk-like 
wind to be consistent with optical and infrared high-resolution 
spectroscopy. This seems to be supported by the present obser- 
vations. In order to be in line with optical blue-shifted emission 
lines and red-shifted absorption lines Oudmaijer suggested the 
scenario of infall of circumstellar material onto the stellar pho- 
tosphere. However, according to Klochkova et al. (1997) the 
concept of accretion does not appear to be unproblematic ei- 
ther. 

Thus, the question which scenario is best suited still ap- 
pears to be a matter of debate. Bispectrum speckle interfer- 
ometry gives important information on the spatial extension of 
the circumstellar shell. It will be in particular the combination 
of different observations - photometry, spectroscopy and high- 
resolution imaging - and their simultaneous modelling, which 
will shed more light on the nature of IRC H- 10 420 that is prob- 
ably being witnessed in its transition to the Wolf-Rayet phase. 
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